Local observables for entanglement witnesses 



G. Mauro D'ArianoQ Chiara MacchiavelloH and Matteo G.A. Pari|] 
Quantum Òptics and Information Group, Istituto Nazionale di Física delia Matèria, Unità di Pavi^ and 
Dipartimento di Fisica "A. Volta", Università di Pavia, Italy 
(Dated: February 1, 2008) 

We present an explícit construction of entanglement witnesses for depolarized states in arbitrary 
fïnite dimension. For infinite dimension we generalize the construction to twin-beams perturbed 
by Gaussian noises in the phase and in the amplitude of the field. We show that entanglement 
detection for all these famílies of states requires only three local measurements. The explícit form 
of the corresponding set of local observables (quorom) needed for entanglement witness is derived. 
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I. INTRODUCTION 

Entanglement plays an essential role in almost all as- 
pects of quantum information theory 0. Entangled 
states are the key ingredients of many quantum pro- 
tocols such as quantum teleportation, quantum dense 
coding, and entanglement-based quantum cryptography. 
However, entanglement can be in general corrupted by 
the interaction with the environment. Therefore, entan- 
gled states that are available for experiments are usu- 
ally mixed states, and it becomes crucial to establish 
whcther or not entanglement has survived the environ- 
mental noise. 

The issue of experimental entanglement detection was 
first addressed for pure states in Ref. Q- More recently, 
in Ref. procedures based on the use of collective mea- 
surements were proposed. Later, in Ref. £| a general 
method to detect entanglement with few local measure- 
ments was presented and optimal schemes were designed 
for two-dimensional systems, bound entangled states and 
entangled states of three qubits. In Ref. [5| a method for 
local detection of nonseparable states has been derived 
for bipartite states in dimension d and to some famílies 
of states of n qubits; it was shown in particular that in 
the bipartite case and for d a prime number the method 
achieves the lower bound of d + 1 measurements derived 
in Ref. Q|. In this paper we extend the approach of 
to depolarised bipartite states in arbitrary dimension, 
and show how entanglement can be efficiently detected 
by identifying the minimal needed set of local observ- 
ables, the so-called quòrum of observables. Moreover, we 
address the problem of entanglement detection for con- 
tinuous variables (CV) and find entanglement witnesses 
(EW) for a twin-beam state (TWB) corrupted by Gaus- 
sian noises, both in the phase and in the amplitude of 
the field. In this case efiïcient homodyne-tomographic 
procedures are analyzed suited to local detection of en- 



tanglement. We found that for all the famílies of states 
that we have considered a rank-four witness operator is 
sufficient to detect entanglement. Notice that this result 
is not in contradiction with the ones derived in Ref. [j| 
because we assume to have more knowledge about the 
family of states. 

The paper is organized as follows. In Sect. [n]we con- 
struct the EW for bipartite depolarized entangled states 
in arbitrary finitc dimension, and give the explícit form 
of the corresponding local quòrum. In Sect. IIIII we ana- 
lyze the case of bipartite CV systems. In particular we 
study the family of twin beam states corrupted by Gaus- 
sian noise, both in the phase and in the amplitude of the 
field, and show how to detect entanglement by employ- 
ing homodyne tomographic techniques. In Sect. I1VI we 
close the paper with a summary of the results and final 
comments. 



II. DEPOLARIZED STATES IN ARBITRARY 
DIMENSION 

In this section we will show how to detect entanglement 
locally for depolarized states in arbitrary finite dimension 
d, namely for the family of states 
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where is any bipartite entangled normalizcd pure 
state of systems with dimension d, I is the dx d identity 
operator and < p < 1. If is a maximally entangled 
state, the states in Eq. Q coincides with the family of 
the socalled isotropic states. 

We will now introduce a more convenient notation. 
Given a bases {\i) <g> \j}} for the Hilbert space Tíi <g> H2 
(with Tíi and H2 generally not isomorphic) , we can write 
any vector I*)) e Hi <g> H2 as 
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The above notation [6j exploits the correspondence be- 
tween states in Hi <8> H2 and Hilbert-Schmidt op- 
erators \& = J2a ^ij\ï)(j\fc° m Hi to H.2- The following 
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relations are an immediate consequence of the definition 
©: 



(3) 
(4) 



where 7? T denotes the transposition of the operator B 
with respect to the chosen basis {\i)}. As mentioned 
above, in the following we will consider only bipartite 
states on Ti <£> H, where Ti has dimension d. 

In this notation the depolarized state QJ takes the 
form 



R= P \nm 
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Let us briefly recali the definition of EW 0, • A state 
p is entangled iff there exists an Hermitian operator W 
such that Tr[Wp} < 0, while Tr[Wp sep } > for all sepa- 
rable states p sep - The operator W is called entanglement 
witness (EW) . For entangled states with non positive par- 
tial transpose (NPT) W can be explicitly constructed as 
W = (|e)(e|) e , where O e denotes the partial transposed 
of O on the second Hilbert space, and |e) is the eigen- 
vector of p e that corresponds to the minimum eigenvalue 
0. Notice that this is not the only method to construct 
entanglement witnesses. Other techniques, working for 
both NPT and PPT entangled states, have been sug- 
gested, as for example, in Refs. |9llïfH|. 

The entangled states of the form have non positive 
partial transpose [ïïl IT^| . Following the approach of £| , 
we will show how to detect entangled states within the 
family JJJ by explicitly deriving EW according to the 
above construction. 

The partial transpose of the state R can be written as 



R e = ® I)E(¥ <g> 7) 
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where E is the swap operator, i.e. E = J2ij |*) 01 ® li) 

As mentioned above, in order to construct a witness 
operator for the family of states Q , we look for the eigen- 
vector of R e corresponding to the minimum eigenvalue. 
Therefore, we can start by writing explicitly the eigen- 
value equation 



R?\A)) = \\A)), 



(7) 



where \A)) is the eigenvector for the eigenvalue A. By 
using the properties J3J and Eq. (JüJ), we can also write 



R e \A))= P \^A T ^*))+c\A)) , 



(8) 



where c = (1 — p)/d 2 , and O* denotes complex conjunc- 
tion of the operator O with respect to the chosen basis 
Therefore, the eigenvalue equation in operatorial 
terms takes the form 



XA = p^>A T ^* + cA 



(9) 



and can be more conveniently written as 

■$A T y*=pA 1 p=(X-c)/p. (10) 

We now use the singular value decomposition of the 
matrix í', namely = XYY\ where X and Y are uni- 
tary operators, while £ is the diagonal operator contain- 
ing the eigenvalues {<Jj} of v$í^ — the so-called singular 
vàlues of ^ — which are conventionally ordered decreas- 
ingly. The above equation then takes the form 



XYY^ A T X*YY T = pA 



(11) 



By multiplying Eq. Ijlll) by X^ on the left and by Y* on 
the right, and upon defining 

B = Y^A T X*, (12) 

Eq. I|ll[) can be written in the compact form 

B T = /i-^SS . (13) 

The last equation can be conveniently expressed by ex- 
plicitly writing its matrix elements as follows 



p l bji(Ji(jj 



By reiterating the above equation one obtains 
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which is fulfillcd for 
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For vàlues of i and j that cannot satisfy Eq. I|16|) we 
necessarily have = 0. We now want to specify the 
form of the operator B corresponding to the minimum 
eigenvalue A. Notice first that for eigenvalues A < c the 
parameter p is negative, and therefore, according to Eq. 
(|14fl . all diagonal elements of B vanish. This is the case in 
particular when the minimum eigenvalue A m is negative. 
We will now explicitly derive the form of B corresponding 
to the minimum eigenvalue A m . Suppose that o\ and 02 
are the two largest elements of E and u\ > 02. Then, 
from Eq. IjlUfl the minimum eigenvalue A m takes the 
form A m = — pa 1O2 + c, and according to Eq. (|15|) the 
matrix elements of the operator B corresponding to A m 
(which we will denote by 7?) are 



-&21 



(17) 



while all the other elements vanish. Therefore, the oper- 
ator B has rank two and takes the explícit form 



B 



1 

-1 


0\ 
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The expression for the operator A corresponding to the 
minimum eigenvalue A m , which we will call A, follows 
from the definition of B in Eq. H12J1 and is given by 

À = XBY T . (19) 

The EW for the family of states QJ can then be derived 

as 

W=(\À)){{À\) 9 . (20) 

Notice that the same form of B is vàlid also for degenerate 
maximum singular value <7i, although in this case the 
solution is not unique. Moreover, an interesting feature 
of the resulting witness operator is that its rank is four, 
independently of the dimension d of the subsystems. We 
also want to point out that the EW for the states Q does 
not depend on the value of p, but only on some a priori 
knowledge about the state namely on the singular 
vàlues of W and on the form of the operators X and Y. 

As an illustration we will consider two explicit exam- 
ples. When \^f)) is a maximally entangled state in di- 
mension d of the form = J^j \ i- e - the opera- 
tor £ is proportional to the identity, with crj = 1/Vd, 
then the operator À corresponding to a state \A)) = 
— \ ji))/y2 can be used to construct a witness oper- 
ator. In this case the state is separable iff p > l/(d + 1). 

As a second example let us consider an initial state 
with Schmidt number two, i.e. ai = a% = 1/V2 and 
(ii = for í > 2. In this case the corresponding EW is 
constructed from \À}) = (|01) - |10))/V2, where |01) and 
1 10) are the basis states related to o\ and 02. The state 
is entangled when p > 2/ (d 2 + 2). 

We will now show how to detect entanglement for the 
family of states (y by measuring only three local observ- 
ables. The matrix A in H19[) can be written as 

A = iX (cr y © 0) Y T , (21) 

where a y is a Pauli matrix (acting between the two levels 
of the two-dimensional subspace spanned by A), de- 
notes the direct sum, and is the null matrix. If P is 
the projection operator over the subspace where A is not 
null, the above expression can be rewritten as 

À = íX'P'Y, y P'Y* , (22) 

where X' = XY T , P' = Y*PY T , and T, y = Y*a y ®0Y T . 
Inserting the above expression in the definition (|20H of 
W = (À ® I)E(Aï <E> I) we have 

W ={X'Y, y ®ï){E 2 ®Q){Y> y X'ï ®I) , (23) 

where Ei is the swap operator for the two-dimensional 
subspace spanned by the support of A. Since one has 

E 2 = \ cr "® cr "' ( 24 ) 

a— t ,x,y,z 

where o t = I, the EW can be finally written as 

W = ~I®I+ \ò a ®o a , (25) 

ct— x : y ,z 



with 

a = X'Yly(T a YlyX'^ . (26) 

As we can see from Eq. (12511 , the witness operator W can 
be measured by performing the measurements of only 
three local observables à a ® o_n, a = x,y,z. This re- 
sult generalizes that of Ref. 4] to arbitrary dimension 
for states of the form QJ: in all cases only three local 
observables are sufficient. 

As mentioned in the introduction, in Ref. :5| a different 
method to detect entanglement of d dimensional states 
has been proposed. This method is vàlid for states of 
the form \ip) = "Y^Z^ a^kk) with > and requires 
the measurement of d + 1 observables. Compared to our 
method, it needs the measurements of a larger number of 
observables, but, on the other hand, it does not require 
the knowledge of the vàlues of the coefncients etfc in the 
density matrix. 



III. PERTURBED TWIN-BEAM IN 
CONTINUOUS VARIABLES 

In this section we address the construction and the 
measurement of EW for CV. At first we have to define 
the families of states we are going to consider. These 
cannot be a trivial generalization of the isotropic states, 
since both maximally entangled states and the identity 
are unphysical states in an infinite dimensional Hilbert 
space. We start from the "maximally" entangled state of 
two CV systems at finite energy, which is given by 

|*)) = * ® I\I}), * = y/1 - \x\ 2 e- xata , \x\ < 1, (27) 

where without loss of generality we will consider x as 
real. Here and in the following, with a\ b\ and a, b we 
will denote the creation and annihilation operators of two 
independent harmònic oscillators, respectively, with com- 
mutations [a, a* ] = [6, 6+] = 1. For the c. m. radiation 
the harmònic oscillators describes two field modes, and 
Eq. (|27|l describes the so-called twin-beam state (TWB) 
obtained by parametric downconversion of the vacuum 
in a nondegenerate optical parametric amplifier. In this 
case n = 2x 2 /(l — x 2 ) represents the average number of 
photons of the TWB. In practice, TWB are the most re- 
liable source of CV entanglement: indeed, experimental 
implementation of quantum information protocols such 
as teleportation, have been obtained using TWB of radi- 
ation. 

Let us now analyze the family of states that are ob- 
tained by perturbing a TWB with a noisy environment. 
We will consider Gaussian noises both in the phase and 
in the amplitude of the field modes. Thermal noise is a 
special case of the present Gaussian displacement noise, 
whereas the noise coming from the addition of a thermal 
state has been considered in |13| . In this case our results 
coincide with the ones given there. 
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The action of a phase-destroying environment on the 
TWB is described by the Master equation 

R = l\2<JaRcJa- {a)a) 2 R- R{cJa) 2 (28) 
+ 2ò t òiïò t ò- (tfbfR- R{tfbf] , (29) 

where R denotes the time derivative of the state R. The 
solution of Eq. for initial condition R = 

can be expressed as 

R(t) = (l-x 2 )Y / x P+q e-^P-^ \pp) (qq\ , (30) 

where we used the abbreviate notation \ij) for |z) ® 
The correlations between the modes are reduced in the 
mixture l|30|) compared to the initial TWB state. How- 
ever, as we will see by explicitly constructing an EW, 
phase-noise never leads to a separable state, i.e. the en- 
tanglemcnt is not destroyed for any valuc of jt. 

In order to obtain an EW for the family R(t) we con- 
struct and diagonalize the partial transpose R e {t) 

R e {t) = (1 - x 2 ) xP+q e-^"^ 2 \pq)(qp\. (31) 

pq 

The eigenvalues equation i? 7 e — is solved by 

A„ = (l-x 2 )x 2n , \ip n ) = \nn) , 
A± m = ±(1 - x 2 )x n+m e- yt( - n - m ') 2 , (32) 

\^nm)) = 7f (\ nm ) ± \ mn )) 

The minimum eigenvalue is given by = —(l — x 2 )xe 1 
corresponding to the eigenvector 

l« = 7f(|01)-|10)). (33) 

The eigenvector IV'qi)) does not depend on 7Í, and thus 
is suitable to build a proper EW for this family of statcs. 
We have 

W = |(hVi>><<V>oil) 9 (34) 
= i(|01)(01| + |10)<10|-|00)(ll|-|ll)<00|). 

The expectation value 

Tr [R(t) W) = Aq X < Ví, x (35) 

is always negative and thus the state R{t) is never sepa- 
rable, for any value of t, and for any value of the initial 
TWB parameter x. In other words, although decreased 
the entanglement is never destroyed by phase-noise. It 
can also be proved 01 that R{t) can be distilled. The re- 
sult in Eq. (IÏÏ5|) proves the conjecture suggested in [Ï5| . 
where the entanglement analysis of a phase-perturbed 
TWB was performed by numerical evaluation of the rel- 
ative entropy of entanglement. 

Let us now consider the family of states obtained by 
perturbing a TWB state by Gaussian amplitude noise, 
namely 

R K = G K ^g K (M(m), (36) 



where for a single mode state p one the map of the Gaus- 
sian noise is given by 

QM = / e-—D(a)ptf(a), (37) 

J 

D(a) = expjaa^ — aa} denoting the displacement op- 
erator. We notice that the operator l|34|l obtained for 
phase-perturbation is an EW also for Gaussian ampli- 
tude noise. Omitting positive factors, we have 

Trliï^^-l + il^^-l + i,. (38) 

Eq. (|38|l says that R K becomes separable if n > 1 — |n _1 , 
a result that can be also obtained by direct check of the 
positivity of the partial transpose (PPT condition) [ït| . 
The family R R , in fact, is composed of Gaussian states, 
for which the PPT condition is necessary and sufhcient 
for separability [Ï7|. It should be mentioned that the 
constructive procedure suggested in Ref. fails to pro- 
vide an EW for the the family R K , in particular it does 
not lead to a state-independent witness. 

In principle, the EW H34|) can be measured by us- 
ing only three observables, as in the finite dimensional 
case. However, there is no feasible implementation of 
the measuring apparatus corresponding to the quòrum 
in the present CV case. Since we are interested only 
in the expectation value of W, we could use quantum 
tomography (for a recent tutorial review on quantum to- 
mography see Ref. ^3)- However, a tomographic deter- 
mination of W is useful only if requires a smaller num- 
ber of observables than those needed for reconstructing 
the full state. Indeed, this is the case for the EW in 
Eq. Q34JI. In fact, for two modes of radiation ai and 
0,2, the expectation value (O) = Tr [RO] of a genèric 
operator O can be obtained by local repeated measure- 
ments of the quadratures Xi^j = ^(aíe 1 ^ 1 -t-etie -1 ^ 1 ) and 
X 24>1 = i(4e 102 +a 2 e" 1 ^) as follows 

(0)=íí ^^(RpKXr^faX^M , (39) 

namely by averaging the over the phases </>i,2 and over 
an ensemble of repeated measurements the function of 
the two quadratures R[0](x±, 4>i; x 2 , 4>2) — so-called es- 
timator or kernel function — depending on the operator 
O. The kernel function for Hilbert-Schmidt operators 
can be obtained directly by means of the trace flsl ] 
R[0]{x 1 ,4> í ;x2 ) <h) = Tr[R(X lcl·l -x 1 )R(X 2cl·2 -X2)0] 
with R(x) = — lim e ^ + ^Reía; + ie)~ 2 . For the operator 
W in Eq. we have 

R[W](xi,<p V ,x 2 ,<f>2) = foo(xi)fn(x2) (40) 
+fn{xi)foo{x 2 ) - 2cos(0i + tf>2)fai(xi)foi(x2) , 

where 

foo(x) = 2$(l,i;-2a; 2 ) 

foi(x) = 4V^$(2,§;-2x 2 ) (41) 

fn(x) = 2[$(l,i;-2a; 2 )-2ci>(2,i;-2 a; 2 )] , 
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and $(a, fe; z) denotes the confluent hypergeometric func- 
tion. Remarkably, R[W] depends only on the sum of the 
two phases 01,2, and shows only a couple of oscillations. 
Therefore, the number of measurements to detect the 
cntanglement witness is much smaller than that needed 
to reconstruct just the first few matrix elements of the 
state, say, in the photon number representation, since the 
number of oscillations of the estimators for such matrix 
elements increases linearly with their photon-number Ín- 
dex. The precision of the tomographic estimation can be 
further improved by adaptive techniques [19). 

If we are allowed to mix the two modes after the per- 
turbation, the characterization of entanglement for the 
family R K can be obtained by measuring a single quadra- 
ture. In fact, for Gaussian states, a necessary and suffi- 
cient condition to have entanglement after a beam split- 
ter is that the two inputs show squeezing (in mutually 
orthogonal directions) [20l |2Ï|. Therefore, if we impinge 
the two modes of the perturbed TWB in a beam splitter, 
and then measure the quadrature X = ~ (a^ + a) on the 
sum mode, we have squeezing if and only if the input 
state is entangled. Therefore, the fluctuation operator 
W = AX 2 - 1/4 = X 2 - (X) 2 - 1/4 is an EW, and its 
expectation value is of course obtained by measuring the 
quadrature X. The analysis is vàlid also when the TWB 
initial parameter x is complex, in which case the phase of 
the quadrature to be measured coincides with the phase 
of x. Obviously, if the mixing of the two modes is not 
possible, one can always reconstruct the above quadra- 
ture locally by quantum tomography. 



IV. CONCLUSIONS 

In this paper we have given an explicit construction 
of EW for depolarised states in arbitrary finite dimen- 
sion. For infinite dimensions, i. e. for CV, we have 
introduced isotropic states as twin-beams perturbed by 
Gaussian noises in the phase or in the amplitude of the 
field, and we have constructed their respective EW as 
well. We have shown that in all cases entanglement de- 
tection needs only a quòrum of three local observables, 
whose explicit form have been derived. For CV it is possi- 
ble to use also homodyne tomography efficiently to detect 
cntanglement, without determining the matrix elements 
of the state. 
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